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ABSTRACT 



The effects of vibrational and chemical nonequilibrium on 
turbine performance were investigated separately. The vibra- 
tional model was taken as a pure nitrogen expansion, and the 
chemical model was taken as the combustion products of a 
stoichiometric mixture of kerosine and air. The loss of per- 
formance of fully frozen flow with respect to the equilibrium 
flow for each model was determined. The extent to which non- 
equilibrium will occur was investigated within limited ranges 
of pressure and temperature. Vibrational nonequilibrium can 
result in losses up to seven percent with respect to equilib- 
rium flow. Vibrational freezing will be virtually complete a 
short distance after the throat of the first stator assembly. 
The losses due to chemical nonequilibrium are insignificant 
compared to the vibrational nonequilibrium losses. 
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I. 



INTRODUCTION 



The aircraft of tomorrow must be capable of operating at 
high Mach number, have large amounts of excess thrust, be 
capable of operating at high altitude, and yet be very eco- 
nomical. As the requirement for these aircraft becomes more 
prevalent, designers of both airframes and power plants must 
strive to match their respective ideas to produce a viable 
design. The airframe specialist is concerned with drag and 
high skin-temperature problems which he hopes to resolve by 
reduction of frontal area and by lighter gross weights. The 
power plant specialist investigates higher turbine inlet tem- 
peratures because of the increased efficiency of the cycle 
with high peak temperatures. The increase in flight Mach 
number induces a rise in the temperature of the cooling air 
available at the turbine. This higher temperature reduces the 
cooling potential of the air. 

The next generations of the turbojet engine will utilize 
nearly double the present turbine inlet temperature by burn- 
ing stoichiometrically . As the turbine inlet temperature is 
increased, the effects of dissociation and vibrational excita- 
tion play a significant role in the energy balance. The jump 
in turbine temperature for stoichiometric combustion will 
produce almost triple the degree of vibrational excitation, a 
large increase in local velocities for the same local Mach 
number within the engine, and a substantial rise in the degree 
of dissociation of the products of combustion. 
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The phenomenon of vibrational energy freezing occurs due 
to the relatively large number of molecular collisions re- 
quired to return the vibrational mode of internal energy to 
equilibrium. Freezing of chemical composition occurs when 
the gas cannot react fast enough to maintain equilibrium with 
a rapidly changing temperature and pressure. Both of these 
phenomena are familiar to the rocket designer. The relatively 
low temperatures and high pressures characteristic to the 
turbojet cycle in the past have not caused these nonequilib- 
rium effects to any significant degree. 

The presence of nonequilibrium in the form of frozen 
vibrational energy and/or frozen chemical composition is very 
probable in the turbine section of an advanced high tempera- 
ture turbojet. The rapid removal of work from the gas in the 
turbine reduces both the temperature and pressure at a rapid 
rate. The vibrational temperature is the independent variable 
in the vibrational energy distribution, and there is a cor- 
responding temperature for both rotation and translation. At 
equilibrium, the rotation, translation, and vibration tempera- 
tures are the same. Comparison of the vibrational, rotational, 
and translational temperatures in the turbine would show a more 
rapid decrease in the latter two relative to the rate of 
vibrational temperature decrease, primarily due to the number 
of collisions required to maintain equilibrium. The magnitude 
of the relaxation time relative to the short transit time in 
the turbine suggests that a significant degree of vibrational 
freezing may result. The effect of lower pressure would be to 
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reduce the collision frequency, thereby increasing the relaxa- 
tion times . 

Evaluation of the loss of available energy due to 
vibrational freezing becomes rather complicated in multi- 
component gases and involves complex-numerical computer pro- 
grams. Stollery and Park [17] and Wilson, et al. [22] give 
examples of these programs. Exact calculations of recombina- 
tion and dissociation rates of a mixture of gases such as 
combustion products also involve complex computer programs. 

The complexity of the entire problem is a result of the num- 
ber of nonlinear equations involved and the lack of an exact 
solution to the flow problem in rotating machines. 

This thesis attempts to apply the nonequilibrium analysis 
of vibrational and chemical freezing to the turbine section 
of a turbojet cycle. The presence and degree of each of these 
problems are investigated, and the effects of changing the 
thermodynamic parameters of the cycle are discussed. 
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II. ANALYSIS 



High, velocity flows frequently create rapid changes in 
the flow variables: temperature, pressure, and density. The 

internal modes of freedom of the gas molecules do not always 
remain In equilibrium or return to equilibrium fast enough to 
keep pace with the rapidly changing macroscopic environment. 
The return to equilibrium of these internal modes is governed 
by internal kinetic equations. The existence of chemical re- 
actions in the flow introduces the chemical kinetic equations. 
There are characteristic times associated with each of these 
processes: relaxation time with the internal kinetic equa- 

tions and reaction times with the chemical kinetic equations. 

Maintenance of equilibrium of chemical reaction and of 
internal degrees of freedom requires a number of molecular 
collisions. It is obvious that true equilibrium flow cannot 
exist. The average time for a unit mass of gas to transit a 
turbine is a prime factor in determining whether or not 
freezing occurs. The conversion of available energy to useful 
work by the turbine is accomplished in a small distance, and 
present turbine technology involves Mach numbers in the 
transonic range. These facts, combined with high local tem- 
peratures, produce a local fluid velocity relative to the 
blade of the order of two to three thousand feet per second. 
The short physical length of the turbine results in an ex- 
tremely low transit time. If the characteristic times of 
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progression of the chemical and molecular phenomena are very 
small compared to the time intervals of interest in the flow, 
the flow may be considered to be in local equilibrium. Con- 
versely, if these characteristic times are much larger than 
the time intervals of interest in the flow, the phenomena 
may be considered frozen. 

Theoretical description of a nonequilibrium flow requires 
incorporation of the kinetic equations describing the processes 
into the usual gas-dynamic conservation equations. Analysis 
and solution of such systems of equations are difficult; in 
all but the simplest problems, numerical methods are required. 
Bray and Appleton [5] and Anderson, et al. [2] , have presented 
detailed analytical solutions to these problems. 

Extensive research into nonequilibrium flows was inspired 
by the advent of space exploration. The flow around re-entry 
vehicles at high altitudes involves complex nonequilibrium 
chemistry [6], In rocket propulsion systems, the problems of 
large expansion ratios and low exit pressures tended to magnify 
the situation. Excellent literature surveys by Read [15], Rich 
and Treanor [16], and Taylor and Bitterman [18] provide a 
myriad of data on vibrational relaxation. Gas dynamic lasers 
employ intentional vibrational nonequilibrium to obtain popula- 
tion inversions [3] . 

The complexity of solution of the expansion of combustion 
products through a turbine should be readily apparent. The 
problem must therefore be simplified in order to obtain 
meaningful conclusions as to the effects of nonequilibrium 
phenomena on the expansion. The interaction of vibrational 
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excitation and chemical reaction is not solved here, but some 
insight into the level of existence of each phenomenon is 
presented. 

The effects and mechanics of vibrational nonequilibrium 
flow are considered separately from the effects and mechanics 
of chemical nonequilibrium flow. 

A. VIBRATIONAL NONEQUILIBRIUM FLOW 

In order to investigate the effects of vibrational non- 
equilibrium in a turbine, a simplified model was utilized. 

The expansion of pure gaseous nitrogen through a turbine was 
considered instead of the more complex case of the actual 
combustion products. Nitrogen was a good model for a number 
of reasons. By using a pure gas, no chemical reaction, other 
than dissociation, was present. The degree of dissociation 
in the temperature range considered was negligible. The prob- 
lem of vibrational freezing could therefore be isolated, and 
the calculation of thermodynamic data for the flow was greatly 
simplified. In addition, the mole fraction of nitrogen is 
large in the gas mixture. 

The following assumptions were expeditious: 

1. Complete (microscopic and macroscopic) equilib- 
rium existed at the exit of the combustion 
chamber. 

2. The gas behaved as an ideal, inviscid fluid. 

3. The turbine expansion was adiabatic. 

4. The flow was considered as quasi one-dimensional. 

5. The velocity of the gas at the combustion chamber 
outlet was considered negligible compared to the 
gas velocity after the turbine exit. 
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A turbine schematic (Figure 1) was defined. At station 
one, corresponding to burner outlet, the nitrogen was con- 
sidered to be in complete equilibrium. Between stations one 
and two, the gas was expanded isentropically through the 
turbine. The work, extracted by the turbine was assumed to be 
drawn strictly from the translational and rotational modes of 
internal energy. These modes were considered to be in equi- 
librium with each other during the entire process due to the 
relatively small number of collisions required for equilibrium. 
This model assumes the independence of the separate modes of 
internal energy. The validity of this procedure may not be 
universal in its application but was utilized here for sim- 
plicity of analysis. The vibrational energy was considered 
frozen in this expansion process. (See Figure 2.) 

Between station two and station three, the gas was allowed 
to return to complete equilibrium in a constant area duct. 

This was accomplished by the redistribution of excess vibra- 
tional energy to the rotational-transational gas through a 
process equivalent to Rayleigh heating. The rotational and 
translational modes were considered to be always in complete 
equilibrium. 

The extraction of work results in a decrease of total 

enthalpy. A hypothetical stagnation temperature of the 

rotational-translational gas at station two, T' , was chosen 

t 2 

to represent a given amount of work extracted by the turbine. 
The prime represents the state of the rotational-translational 
gas. The ratio of specific heats, y, was taken as a constant 
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and equal to that obtained by classical equipartition con- 
cepts for a diatomic gas with only rotational and transla- 
tional modes. This value is y = 1.4. The selection of a 
burner outlet stagnation pressure and temperature and a Mach 
number at station two, based on the chosen rotational - 
translational (RT) gas temperature, enabled the determination 
of the station two stagnation pressure through the use of the 
ideal gas isentropic relationship in equation (1) . 




Use of an isentropic procedure is consistent with the desire 
to isolate vibrational nonequilibrium. 

Downstream of station two, the gas was subjected to a 
Rayleigh heating process due to the relaxation of the vibra- 
tiona 1 mode . 



Q + C D Ti 
r P RT 2 





(2) 



Utilizing the definition of stagnation enthalpy, 



<2 r + H 2 = H 3 (3) 

where Q represents the energy added to the RT gas from the 
vibrational mode. 

Equation (3) may be rewritten as: 



Q„ = h* - hi = Cp (t; - T/ ) (4) 

RT 2 

The ideal gas isentropic relationship was used to find static 
temperature, T^, which was needed in the iteration procedure. 
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( 5 ) 




The determination of the Mach number at point three was then 
required. This value was determined by an iterative process 
matching the Rayleigh heat added to the classic definition of 
enthalpy. The enthalpy was considered as 



h (T) 




C p dT 



o 

Since the internal energy modes were considered separable 
specific heat may be expressed as 



(6) 

the 



Cp Cp + Cp + Cp 

R T V 



= C. 



+ C r 



(7) 



RT V 

Using equation (7) , equation (6) may be expressed as 



h (T) 



■f 



C dT + 
P RT 



I 



C p dT 
v 



(8) 



o o 

The energy added to the equilibrium RT gas in the Rayleigh 
heating model must be identically equal to the difference in 
the vibrational energy of the complete gas between stations 
one and three. This difference was taken thusly 



h(T 1 ) 



r T i r T i 

= J J ' 



Cp dT + Cp dT 

‘v J ‘RT 

o o 



(9) 



Therefore 



h(T 1 ) 



■ Jl\ 



dT + C-, (T.-T ) 

P RT 1 ° 
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Similarly 



h(T 3 ) 



■jj 



C P dT + C P ot < T 3- T o> 
v RT 



( 10 ) 



Taking the difference between equations (9) and (10) 



hCT x ) - h(T 3 ) 



■ r 



Cp dT + C (T,-T-) 
v RT 



(ID 



Now define, Qp, as the difference in the vibrational energy 
between stations two and three. 

• T, 



■ i 



C p dT 
v 



( 12 ) 



Substituting equations (11) into equation (12) 



Q 



D 



h(T ) - h(T ) - C (T.-T ) 

r RT 



(13) 



The iteration was satisfied when Q^_ of equation (4) was equal 

to of equation (12) . The actual solution for a given set 

of input parameters was obtained by entering gas tables [11] 

with the chosen M ' . With the already chosen value of T' , the 

Z t 2 



values of T* and T* were obtained. A Mach number at station 



3' was assumed, and equation (4) was solved in the following 
form. 



= T* C 
1 t 



RT 



T i. 



T t. 



T* 



T* 



(14) 



Utilizing information from gas tables [11] the temperature, 
T^, was also obtained and equation (12) was evaluated. This 
process was repeated until Q r and Q D were within 0.5% of each 
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other. Satisfaction of the iteration process allowed calcula- 
tion of stagnation temperature and pressure at station three. 

The investigation of performance lost due to vibrational 
freezing required a reference quantity. This reference quan- 
tity was taken as the turbine work produced by the expansion 
of the gas with full equilibrium between the modes of vibra- 
tion, rotation, and translation. For both equilibrium and 
nonequilibrium cases, the turbine pressure ratios were assumed 
identical. 

Efficiency was defined as the ratio of actual work to 
ideal work where the actual work was considered as the work 
obtained from the vibrationally frozen model and the ideal 
work as that obtained from an equilibrium gas expanding to the 
same stagnation pressure. 



where T* is the ideal stagnation temperature at the stagnation 
pressure of station three. 

The stagnation temperature for an ideal turbine evaluated 
.at the stagnation pressure of station three was determined by 
isentropic relationships. By definition. 



n 



vibrationally frozen work 
v equilibrium work 



( 15 ) 



H (T ) - H (T ) 
r l C 3 



n v H (T ) - H("T* ) 

t x t 3 



( 16 ) 



S(T. , P. ) 
r l C 1 




( 17 ) 



and 



S(T, P) 



S° (T) - R £n (P) 



( 18 ) 
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where S° (T) is the entropy of the gas at a given temperature 
and a pressure of one atmosphere. Using equation (18) in 
equation (17) 

S° (T ) - R S-n P = S° (T* ) - R Jin P* (19) 

r l r l c 3 t 3 

It was chosen that P* = P' and the values of T and P 

r 3 *3 r l fc l 



were assumed. Therefore, an equation for the entropy at the 
desired temperature was obtained. 



S° (T* ) = S°(T. ) + R In 

fc 3 fc l 

Solution of this equation yielded T* 

3 



P £. 



( 20 ) 



which was utilized to 



compute values of efficiency from equation (16) . This effi- 
ciency represents a loss of performance resulting from non- 
equilibrium flow. 

A method to identify the significance of vibrational non- 
equilibrium had been established, and the question of actual 
presence of this phenomenon in a high temperature turbine 
remained. 

The complexity of the problem of vibrational relaxation 
precludes actual solution of the theoretical equations within 
the scope of this thesis, and experimental results of other 
researchers will be utilized in the ultimate discussion of 
the effects of this phenomenon on turbine flow. For complete 
understanding it is prudent, however, to develop the princi- 
ples involved. 

The nitrogen gas was considered to be a system of harmonic 
oScLllators . This system of diatomic molecules was assumed to 
be capable of collisional interchange of energy with the 
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translational and rotational degrees of freedom of a heat 
bath at temperature T. The heat bath was provided by the 
translational and rotational degrees of freedom of the pure 
gas. At equilibrium/ the distribution of oscillators over 
the energy states was assumed to be given by the Boltzmann 
distribution . 

The majority of data on vibrational nonequilibrium, both 
theoretical and experimental, has been developed for the lag 
of vibrational equilibrium when vibrational degrees of free- 
dom are being excited. Rich and Treanor [16] point out that 
there is evidence that the relaxation time for de-exitation 
may be one-fifth or smaller relative to the exitation relaxa- 
tion time. The relaxation time is a measure of the rate of 
return to equilibrium after a disturbance. 

The equation for the time rate of change of the number of 
molecules in a particular energy state may be developed using 
statistical thermodynamics and quantum mechanics. The rate 
of transition from one state to another is dependent upon the 
number of collisions per unit time between the oscillators 
and the molecules of the heat bath that possess sufficient 
energy to cause transition and on the probability that a 
single collision will result in transition. 

For gas dynamic purposes the total vibrational energy is 
more useful than the molecular energy level population in 
each energy level. This total energy is a summation over all 
possible levels of the product of the energy in a particular 
level and the number of molecules in that level. Manipulation 
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of the rate of change of the molecular level population equa- 
tion produces an equation for the rate of change of total 
vibrational energy. Vincenti and Kruger [20] express this as 



dE 



v 



E* - E 
v v 



dt 



( 21 ) 



where E v is total vibrational energy, E* is total vibrational 
energy at equilibrium, and x is the relaxation time. 

This equation shows that the vibrational energy of the 
system will always move toward equilibrium and that the rate 
of this movement at any instant is linearly related to the 
degree of departure from equilibrium at that instant. Even 
though this equation assumed a Boltzmann distribution only at 
equilibrium and no restriction is placed on the amount of 
departure of vibrational energy from equilibrium, Vincenti and 
Kruger [20] limit the equation to small departures from equi- 
librium. The limitation to small departures is a result of 
the anharmonicity of actual molecules. 

The relaxation time, t , is a function of the temperature 
and pressure of the heat bath and of the characteristics of 
the vibrating species. For practical purposes it was assumed 
that equation (21) was valid for any number of excited mole- 
cules and could be applied to a time varying heat bath as long 

as E* and x were taken as functions of the instantaneous state 

v 

of the heat bath. Under these assumptions, x should be tech- 
nically referenced as a local relaxation time where numerical 
integration could produce a solution for a particular time 
span . 
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For sufficiently low temperatures the original Landau and 



Teller Theory [20] may be approximated by: 



C , K 2x 
p exp (^-) 



1/3 



T 



( 22 ) 



where K is a constant dependent upon the physical properties 
of the species, and C is an empirical constant dependent upon 
the type of heat bath molecules. 

In the present theoretical state, the evaluation of t must 
be accomplished by experimental measurements interpreted in 
the light of the foregoing theoretical results. The methods 
of ultrasonic absorption, shock tube interferometry, spectro- 
scopic data, and pressure measurement, have been used experi- 
mentally to produce the values of the constants for use in 
equation (22). These values are presented in Table II. 

Taylor and Bitterman [18] present numerous graphs of both 
t and the probability per collision of an exchange of energy 
between the translational and vibrational modes of a molecule. 

Anderson, et al. [2] , presents correlations between re- 
laxation times, temperature and pressure for nitrogen and 
various collision partners. These correlations were utilized 
extensively and are therefore presented here. 




93 (T 1/3 ) - 4.61 



(23) 



log io ^ tp) n 2 -co 2 




(24) 



l°g 10 (TP) 




15.4 (T 1/3 ) + 0.722 



(25) 
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Bray [4] suggests that once an appreciable deviation from 
equilibrium has occurred that return to equilibrium in a 
nozzle is unlikely. Therefore, the investigation was cen- 
tered on the stator nozzle block of the turbine. A stator 
passage with area change from throat to exit given by equation 
(26) and shown in Figure 3 was chosen, 

A/A* = 0.15x - 0.0 lx 2 (26) 

where A* was the area at the throat and x was the distance 
along the mid-span streamline measured from the throat. While 
the choice of this particular formula for change of passage 
area was arbitrary and not particularly critical, this choice 
is representative of modern turbines. 

Knowledge of the area change and selection of a represen- 
tative y enabled determination of Mach number, and thermo- 
dynamic data from gas tables as presented in Table III. 
Selection of a total pressure and temperature at the burner 
outlet permitted evaluation of vibrational relaxation times 
at the throat of the stator with area change given by equation 
(26) for nitrogen to nitrogen collisions using equations (22) , 
(23), or the graphs of Taylor and Bitterman [18], (See Table 
IV.) 

As stated previously, the calculated relaxation times had 
to be compared to a gas transit time. The gas transit (dwell) 
times were determined by dividing a characteristic length by 
the appropriate velocity. Under this concept, a total transit 
time was obtained by dividing the distance from point A to 
point B along the mid-span streamline by the mean velocity 
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between those two points. A local dwell time at a data point 
was obtained by dividing the distance between data points by 
the velocity at the upstream point. 

Numerical evaluation of equation (21) was performed 
utilizing calculated relaxation times and the method of Hurle , 
Russo, and Hall [9]. The method is illustrated in Appendix C. 
This evaluation illustrated the actual amount of vibrational 
energy as compared to the equilibrium amount. 

An aircraft gas turbine using hydrocarbon fuels will have 
many molecular species in the combustion products. It is 
known that some of these species are more effective for vibra- 
tional de-excitation. Consequently, a composite relaxation 
time involving several molecules was calculated. The molecular 
species in the combustion products were screened to ascertain 
the molecules most effective for vibrational de-excitation. 

This screening indicated ^0 and CC^ should be included. A 
relaxation time was thus calculated for nitrogen with three 
different collision partners, ^ , CC> 2 , us i- n 9 equations 

(23) , (24) , and (25) respectively (Table V) . The total rate 

of change of vibrational energy of the nitrogen was considered 
to be the sum of three distinct quantities corresponding to 
the change induced by each of the collision partners. 

dE dE dE dE 

-rr 1 = A + B + C (27) 

a total ^2 _N 2 ^ 1 2“ C0 2 ^2 _H 2° 

Each of the terms on the right-hand side of this equation was 

expressed in the form of equation (21) and substituted to form 

equation (.2 8) 
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dE 



v 



dt. 



E *-E E*-E 

v v v v 

= A — — + B — — 



+ C 



E *-E 
v v 



(28) 



n 2 -n 2 



n 2 -co 2 



n 2 -h 2 o 



'total 

The coefficients (A, B, C) were taken to be the mole fractions 
of the collision partner. Since the relaxation times with N 2 
and CC> 2 are equal, equation (2 8) may be expressed as 

dE 



v 



dt 



= (E*-E ) 

' v v' 



1 - 



v + V 



N 2~ N 2 



n 2 -h 2 o 



(29) 



where ^ is the mole fraction of water in a gas mixture. 
Comparing equation (29) with equation (21) 



1 

T 

composite 





t n 2 -h 2 o 



(30a) 



The extension of equation (30a) to systems with a larger num- 
ber of components is obvious. Note that E v is the vibrational 
energy of nitrogen only. Values of x composite are in Table 

VI. 

Consider a large number of N*, where N* denotes vibra- 
tionally excited N 2 . If the number of N* exceeds the equilib- 
rium value, there will be a decrease in the number of N* as 
time progresses. The N* may be de-excited by collisions with 
H 2 0 or C0 2 or N 2 . The collision frequency with H 2 0, CC> 2 , or 
N 2 depends upon the mole fraction. The probability that a 
given collision will result in de-excitation depends on the 
partners; i.e. on C and K 2 in equation (22). Using this 
physical picture, equations (27), (28), and (29) follow. In 

equation (29), each t is evaluated at a pressure equal to that 
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exerted by the gas mixture. Alternately, an expression for 
x composite would be 



1 



1 



1 



(30b) 



T 

composite 



x 




x 




where partial pressures are used to evaluate x for each 
species. The analogy of £he equation for x composite to the 
formula for resistances in parallel is apparent. 

B. CHEMICAL NONEQUILIBRIUM FLOW 

In' order to investigate the effects and mechanics of 
chemical nonequilibrium on the flow in a turbine, the products 
of combustion of a stoichiometric reaction between air and a 
kerosine-type hydrocarbon fuel were considered and a turbine 
expansion calculation performed. The choice of a stoichio- 
metric mixture ratio was indicated by the current trends and 
research goals in the turbojet industry. The internal degrees 
of freedom of all components of the gas were assumed to be in 
equilibrium throughout the expansion. Even though there is a 
known relationship between vibrational modes and chemical re- 
action, this assumption was necessary for analytical purposes. 

The gas was taken as an inviscid, diffusionless fluid with 
quasi one -dimensional expansion through an ideal turbine. 

The process was considered to be adiabatic. The gas was ex- 
panded first as chemically frozen and then in full chemical 
equilibrium. The stations in the process were as illustrated 
in Figure 8. The concepts used in this section closely follow 
Penner's [13] development for chemically frozen flow.. 
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In chemically frozen flow, the average molecular weight 
obviously remains constant; therefore, the work derived from 
this expansion can be expressed as 



Work 



frozen 



n 

T X. (H. (T. 

i-i 1 1 4 



i - 



V T t >> 



(31) 



where is the molar fraction of species i at station 1, n is 

the number of .-.species , present , and H^(T) is the total molar 

enthalpy of species i at temperature T. 

Given values of P , T , and T , the total pressure at 

r l t l C 2 



station two was determined from the isentropic relationship 



= P. 



Y* 

y*-l 



(32) 



where y* is an appropriately defined average specific heat 
ratio during expansion. Penner [13] gives this ratio as: 



Y* 



l X. (S° (T. ) 

i=l 1 r l 



S°(T. )) 

2 



n _ iti 

l X (S° (T ) - S° (T ) ) - R Jin =-^= 
i=l 1 r l r 2 X t 2 



(33) 



The chemical nonequilibrium efficiency, n / was defined 

c 

as the ratio of the work extracted from chemically frozen gas 
to the work extracted from equilibrium gas expanded to the 
same total pressure. This efficiency is 



H 1 ~ H 2 
Mi 

H 1 " H 3 ~ 
M 3 



(34) 
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where 





n 






H 1 = 


i x i< T t - 

i=l 1 C 1 


P. ) H . (T. ) 
C 1 1 C 1 


(35) 




n 






II 

<N 


l x i( T t - 

i=l 1 Z 1 


P. ) H . (T. ) 

r l 1 Z 2 


(36) 




n 






H 3 " 


I X i(T t » 
i=l 1 r 3 


P+ )H. (T. ) 

3 Z 3 


(37) 




n 






II 


l X i(T t . 
i=l ^ 


P tl )m i 


(38) 




n 






M = 

3 


l X^T , 
i=l 1 C 3 


p t 3 )ra i 


(39) 



Once the quantities in the numerator of equation (34) 
were established as functions of known or given quantities, 
attention was focused on the denominator. 

For equilibrium flow, the analysis was restricted to a 
given weight of gas mixture due to the temperature dependence 
of the average molecular weight of the mixture. For purposes 
of comparison it was convenient to choose the initial mole- 
cular weight of the gas mixtures as the characteristic 
weight; thus the ideal or equilibrium flow work was: 



Work 



equi 



where 




(40) 



H 



1 



l X. (T, , Pt,) [aH° + H. (T. ) 

i=l X Z 1 1 L i 1 Z 1 



E. (T 

X 



o'] 



(41) 
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and 



H. 



I X. (T. , P. ) [AH 0 + H. (T ) - H. (T j| (42) 

i=l 1 c 3 3 L i 1 r l 1 °J 



At this point, all the necessary quantities for evaluation of 

H were known except T which could be determined from 
C t 3 

isentropic considerations. 

.Isentropic flow occurred for the total weight of the gas 
mixture. Since was chosen as the fixed weight for analysis, 
the isentropic relationship was: 

M, 



1 7T, 



S(T t , P t ) = -±- S(T t , P fc ) 



(43) 



M. 



S (T,P) is the molar entropy of an ideal gas mixture at tem- 
perature T and pressure P and may be expressed as 
n _ n 

S (T , P ) =' l X . (T , P) S ? (T) - RUn P + £ X. (T,P) £nX. (T ,P) ) 

i=l 11 i=l 1 1 

(44) 

where S° (T) is the molar entropy at temperature T and one 
atmosphere of pressure. Equations (44) and (44) were combined 
to yield 
n 



X. (T , P )fs?(T ) - R £n P - R S,n X. (T , P )1 

i=l ^1 r l L 1 1 1 1 J 

M. n p _ 

= zr l X. (T. ,P. ) s? (T. ) - R JJ-n P 

M- i=l 3 3 L 3 t 3 



- R Jin X. (T , P ) 
x t 3 t 3 



(45) 



This equation was used to solve for T by iteration. Once 

fc 3 

this has been accomplished, evaluation of equation (34) was 
possible . 
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The mechanics of the chemically reacting flow were then 
investigated. In the combustion chamber of a turbojet, a 
significant portion of the chemical energy released is used 
to break chemical bonds producing more elementary species. 

This energy is unavailable to do work unless it is first re- 
leased through a recombination of the dissociated combustion 
products. The rate of each chemical reaction must be compared 
to a dwell time of the flow at any particular stage in order 
to determine the extent of the unavailability of energy. 

Consider the typical three body reaction 

r f 

A + B + M+ AB + M (46) 

r b 

where M is a third body which removes energy during the re- 
combination and r - and r, are reaction rates in the forward 

t b 

(recombination) and backward (dissociation) directions 
respectively. This equation represents two separate reactions 
with the connection that r^ must equal r^ at equilibrium. As 
the temperature and pressure move toward lower values during 
the expansion process, the reaction in equation (46) should 
move in the forward direction. If the reaction rate is not 
fast enough to maintain equilibrium conditions, recombination 
will lag, and the gas driving the turbine will have dissociated 
products whose dissociation energy has not contributed to the 
kinetic energy of the flow. The lag of recombination results 
in a reduction of performance. 

Equilibrium flow calculations involve the simultaneous 
solution of the conservation equations for mass, energy, 
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momentum, and species, plus the equation of state. This sys- 
tem is simplified for frozen flow because the species composi- 
tion is constant. Computer programs for solution of 
equilibrium and nonequilibrium flows are presented by Gordon 
and Zeleznik [8] and Zupnik, et al., [23], 

For relaxing flow in a turbine the problem is complex 
because the nonlinear kinetic differential equations of 
species concentration and the conservation of atom equations 
are now involved. Solution of the chemical relaxation problem 
is obviously desirable but beyond the scope of the present in- 
vestigation. Penner [12] and Bray [4] have proposed criteria 
to determine when the shift may be made from equilibrium to 
nonequilibrium analysis. Application of these criteria by 
various researchers has produced varied degrees of success. 

Newhall [10] presents twenty -eight elementary reactions 
(Table X) that may be considered significant in the flow of 
combustion products resulting from a mixture of a hydrocarbon 
fuel and air. These equations may be integrated numerically 
with the coupled energy equation to produce concentration- 
time histories for each species during expansion, but this 
was beyond the scope of this thesis. 

The science of chemical kinetics, the study of the rates 
of chemical reactions, is complex and difficult, both experi- 
mentally and theoretically. A discussion of both simple col- 
lisional theory and transition state theory is given in Pratt 
[14] and emphasizes the heavy reliance, except in the simplest 
reactions, upon empirical data. 
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The chemical reaction rate may be considered to be pro- 
portional to the total number of collisions between reacting 
chemical species and to the fraction of this total in which 
sufficient energy is available for the reaction to occur. In 
flow systems the rate of reaction is also dependent on the 
mass and energy transport mechanisms. The collision rate is 
dependent upon both concentration and effective collision 
cross section of the molecules involved. 

The most widely utilized approach to the reaction rate 
problem is through the use of the Arrhenius equation 



K = A (T) exp ( -E /RT) (47) 

cl 

where A is the pre-exponential factor and E is the experi- 

a 

mental activation energy. The values of both A and E are 
obtained from experimental plots of £n K vs l/T. Both A and 
E may be functions of temperature. 

For elementary reactions, such as those that make up the 
complex reactions in the combustion products, the rate con-_ 
stant, K, from the Arrhenius equation (47) may be employed as 
follows: The rate of change in the concentration of a product 

AB produced by a reaction such as equation (46) . 



A + B + MjAB + M 



(46) 



is equal to 



d (AB) 
dt 



K (A) (B) (M) 



(48) 



for an elementary reaction where the parentheses signify 
concentration in moles per cc. 
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Each of the equations in Table X had corresponding for- 
ward and backward rate constants from Table XI. Each of 
these rate constants was calculated at the throat of the 
turbine stator model introduced previously (Table XII) . The 
mole fractions of the species in the combustion products were 
calculated at the stator throat using Banes, et al. , [3] and 

are presented in Tables XIII and XIV. A conversion factor 
that enabled the use of mole fractions instead of molar con- 
centrations in equation (48) was developed and is presented in 
Appendix D. The reaction rates for the equations considered 
were calculated at the stator throat and are tabulated in 
Table XV. The total reaction rate for each specie was deter- 
mined by a summation of the rate of destruction or creation of 
the specie from the applicable equation. If the particular 
specie appeared on the left-hand side of the reaction equa- 
tion in Table X it was considered as being destroyed at the 
forward rate and produced at the backward. In Table XI, the 
total reaction rate at the stator throat for each specie is 
presented. 

The determination of the presence of chemically frozen 
flow in the turbine required establishment of a characteris- 
tic reaction time. The total reaction rate, R^ , of a specie 

-1 -3 

r, had units of moles sec cm . Multiplying by the 

3 

specie molecular weight produced units of grams of i per cm 
per sec, which yields a rate of change of density. 

(49) 



R. m. 
r r 



dp i 

dt 



By taking the gas average density at the throat, p^, and 
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dividing both sides of equation (49) by p\ , a rate of change 
of mass fraction is obtained for small mixture density 
changes . 



R. m. 

1 x 



1_ 

- dt 
P < 



dp. dY. 

1 „ 1 

“ dt 



By definition 



m. 

-i- X. 
M 1 



For sufficiently small changes 



, v Y. - Y. 

^i ~ _^2 H 

dt t 2 - t x 



(50) 



(51) 



(52) 



Letting t 2 - t^ = t equation (50) may be rewritten as 



or 



R. m. 
x x 



dY. 

x 

dt 



Y. - Y. 
X 2 X 1 

fc R. 

1 



'R. 

1 



m. m. 

X. ^ - X. ^ 

1 2 M 2 X 1 M x 

r! in! 

x x 



(53) 



Rearranging 



'R. 

x 



R. 



X. 

x. 



M, 



X. 

x. 



M, 



(54) 



‘2 “1 

Evaluation of equation (54) provided a basis for comparing 
reaction times with dwell times. 
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III. RESULTS AND DISCUSSION 



Values of the vibrational efficiency defined by equation 
(16) were calculated using stagnation conditions of 2200 °K 
and 20 Atmospheres. These values are presented in Table I. 

The data in Table I are plotted in Figure 4 as loss percentage 
against the change of total temperature of the RT gas from 
stations 1 and 2. See Figure 1 for station numbering. Selec- 
tion of this abscissa illustrates the variance of the vibra- 
tional efficiency with the amount of work extracted. 

The data for different Mach numbers indicates that this 
efficiency as defined is a definite function of turbine exit 
Mach number. It is known that for a given static temperature 
there is a known amount of energy stored in the vibrational 
mode of internal energy at equilibrium. As the exit Mach num- 
ber increases, the kinetic energy extracts a larger portion of 
the rotational and translational energy depicted in Figure 2. 
Therefore, the static temperature of the RT gas is lower and 
the amount of equilibrium vibrational energy is lower. The 
total vibrational energy at the turbine exit is the same as 
the equilibrium vibrational energy at the burner outlet. As 
the Mach number increases, the RT gas static temperature de- 
creases and the fraction of the total vibrational energy at 
the exit that corresponds to the RT gas temperature decreases. 
Therefore, the flow experiences a larger departure from equi- 
librium and an attendant increase in the losses. 
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The apparent minima of curves in Figure 4 as indicated by 
the dashed line is a result of the mechanism of work removal. 
The work was extracted from a RT gas which had both modes 
fully excited in the temperature range considered. This full 
excitation made the work extraction a linear function of tem- 
perature. The equilibrium work in the denominator of equation 
(15) was drawn from all modes of the gas. At the temperatures 
considered, vibrational specific heat is nonlinear and is the 
apparent cause of the minima. 

There were three sources of vibrational relaxation times: 
equation (22) , the graphs of Taylor and Bitterman [18] , and 
the equations of Anderson, et al., [2]. In Table IV, relaxa- 
tion times for nitrogen-nitrogen collisions are calculated by 
the three methods. The values are for an area change as de- 
fined in equation (25) and stagnation conditions of 2200°K and 
20 atmospheres. The variation in the values of relaxation 
time when calculated by the three methods indicates that the 
data basis for each calculating method is different. The 
graphs of Taylor and Bitterman [18] had insufficient range 
.for this analysis. The values obtained through the use of 
equation (22) seem unusually small and therefore were not used. 
The Anderson, et al. , [2] data is a later development with 

information for many different collision partners, and its 
use allowed more flexibility of investigation. 

The relaxation times for nitrogen as a function of stator 
area ratio are given in Table V for a stagnation temperature 
of 2200 °K and various stagnation pressure ranging from 5 to 
20 atmospheres. Equation (26) was used to calculate area ratio 
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as a function of distance along the midspan streamline. These 
data are plotted in Figure 5 with the area ratio as the 
abscissa. The area ratio values are of significance only in 
that they define data collection points. The pressure de- 
pendence of the relaxation is aptly demonstrated in this 
figure. The lower pressure results in a lower collision fre- 
quency and thereby slows the relaxation. The curve of relaxa- 
tion time for nitrogen-water collisions emphasizes the 
catalytic effect of water on the relaxation process. The 
presence of water in the combustion products of a hydrocarbon 
and air mixture helps to delay the onset of vibrational 
freezing. 

The dwell time as plotted in Figure 5 represents the re- 
ciprocal of the local velocity in units of seconds per inch 
multiplied by a factor of ten. The order of magnitude of the 
dwell time was raised because exact numerical solution of 
equation (21) indicated vibrational freezing is significant 
at A/A* = 1.14. The factor of ten multiplication results in 
the dwell time and relaxation time curve intersection agree- 
ing with that solution. The multiplication factor makes the 
dwell time longer and more pessimistic. The intersection of 
the dwell time curve with the relaxation time curves represents 
a position downstream of which the flow may be considered to 
be vibrationally frozen. For the pure nitrogen, this point 
occurs at a low area ratio indicating that the flow is prob- 
ably fully frozen for practical calculations. Once the flow 
has frozen, Bray [4] indicates that it will remain frozen or 
at least in some form of nonequilibrium. The fact that 
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freezing is indicated early in the stator section of the 
turbine implies that the flow should be considered as vibra- 
tionally frozen throughout the turbine performance 
calculations . 

Equation (21) was evaluated for various conditions and the 
data plotted in Figures 6 and 7. Both of these figures dra- 
matically demonstrate the departure from vibrational equilib- 
rium that occurs at a low area ratio. The effects of other 
collision partners are shown in Figure 7 using data from 
Tables VII and VIII. The upward movement of the curve in 
Figure 7 when the midspan streamline length is reduced is a 
result of a lower gas transit time. 

Values of the chemical nonequilibrium efficiency were com- 
puted from equation (34) for stagnation conditions of 2200°K 
and 20 atmospheres of pressure at the burner outlet. These 
values are presented in Table IX. Data from Table IX were 
plotted as loss percentage in Figure 9 with the change of 
total temperature from burner outlet to turbine exit as the 
abscissa. Here again the total temperature change is an indi- 
cator of the amount of work extracted. The losses due to 
chemical freezing in this stoichiometric-mixture expansion are 
rather low. For the temperature range considered, the amount 
of nitrogen dissociation is negligible. This is significant 
because nitrogen comprises nearly three-quarters of the gas 
mixture. For the stoichiometric mixture there is very little 
oxygen remaining to dissociate after the reaction. Consider- 
ing the size of the losses, the tendencies of the curve are 
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subject to question. The precision of the thermochemical 
data utilized for the calculation would have significant ef- 
fect on the results. 

The characteristic reaction time for each specie was 
calculated from equation (54) using the data for the stator 
throat and stator exit at total temperature of 2200 °K and 
total pressures of 15 and 5 atmospheres. These times are in 
Table XVIII. Only the reactions considered to be the fastest 
were utilized in the analysis. The influence of slower reac- 
tions upon the fast reactions was not analyzed. Each of the 
times was analyzed and a classification was performed. The 
gas transit time from stator throat to exit was compared to 
each specie reaction time, and the specie was placed in an 
appropriate section of Table XIX. 

A negative reaction time for a specie indicated that the 
predominate or fast reactions involving that specie were away 
from equilibrium for the particular conditions of temperature, 
pressure and concentrations. The real life case, including 
reactions not considered here, would, for long time periods, 
cause the reaction times to become positive. The negative 
reaction time species were classified as frozen. A long or 
short positive reaction time compared to the transit time 
classified the specie as frozen or equilibrium respectively. 
Specie reaction times of the same order of magnitude as the 
transit times were classified as unknown. The approximation 
used to calculate the reaction times was not sufficiently 
accurate to make a definite statement as to the freezing of 
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the particular specie. Table XIX definitely indicates the 
presence of chemically near-frozen flow. 

Several aspects of the current state of the art in turbo- 
jet development tend to reinforce the gravity of the losses 
calculated in this paper. As the aircraft operate at higher 
and higher altitudes, the total pressure at the burner outlet 
will decrease. This will insure the presence of vibrational 
and chemical nonequilibrium as illustrated in Figure 5, 6, 7, 
and Table XIX. The use of hydrogen fuel would produce large 
amounts of water thereby delaying the vibrational freezing 
effects. The magnitude of the burner outlet temperature will 
vary as the ability of the turbine blade cooling method to 
handle high inlet temperature, and this variation will effect 
the degree of nonequilibrium present. 
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IV. CONCLUSIONS 



The vibrational freezing of the flow in the turbine 
section of an advanced turbojet using a stoichiometric mix- 
ture of kerosine and air can result in losses of up to seven 
percent when compared to equilibrium flow. 

The vibrational freezing of the flow will be virtually 
complete a short distance after the throat of the first 
stator assembly. 

The losses due to chemical freezing are insignificant 
compared to the vibrational freezing losses. 
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APPENDIX A 



TABLE I 

VIBRATIONAL EFFICIENCY WITH TOTAL TEMPERATURE OF 2200°K 
AND TOTAL PRESSURE OF 20 ATM 



T^, °K 


M 2 “ - 1 




M 2 ' -3 




M 2 = * 5 


1500 




96.37 




95.21 




92.66 


1550 




96.85 




95. 40 




92.63 


1600 




97.00 




95.52 




92.49 


1650 




97.27 




95.95 




92.43 


1700 




97.71 




95.89 




92.02 


1750 




97.90 




95.80 




91.62 


1800 




97.98 




95.80 




91.03 


1850 




97.81 




95.43 




90.51 


1900 




98.14 




95.14 




89.46 






TABLE II 








VALUES OF CONSTANTS 


FOR USE IN EQUATION (22) 






FROM 


VINCENTI 


AND KRUGER [20] 






Species 


Heat Bath 
| Molecules 


C 

| Atm-microsec | 


K 2 ° K | 


Approximate Range 
of T °K 


°2 


°2 


5.42 x 


10 " 5 


2.95 x 10 6 


800 


- 3200 


°2 


Ar 


3.5 x 


10" 4 


2.95 x 10 6 


1300 


- 4300 


N 2 


N 2 


7.12 x 


10 " 3 


1.91 x 10 6 


800 


- 6000 


NO 


NO 


4.86 x 


10" 3 


1.37 x 10 5 


1500 


- 3000 


NO 


NO 


6.16 x 


10- 1 


1.37 x 10 5 


1500 


- 4600 
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TABLE III 



STATOR AREA CHANGE WITH ASSOCIATED IDEAL GAS DATA 







FROM WANG, ET AL. 
y = 1.26 


[21] 




A/A* 


M 


T/T t 


p/p t 


Velocity , 
M/sec 


Dwell Time, 
sec/in 


-.1 


1.0 


. 885 


.5531 


854 


29.7 


1.14 


1.43 


.791 


.3191 


115 4 


22 


1.27 


1.59 


.7526 


.2523 


1252 


20.3 


1.37 


1.69 


. 7292 


.217 


1310 


19.4 


1.46 


1.77 


.7106 


.1925 


1354 


18.8 


1.51 


1.81 


. 702 


. 18 


13 76 


18.5 


1.55 


1.83 


.696 


.172 


1389 


18.3 



TABLE IV 

COMPARISON BETWEEN VARIOUS METHODS OF OBTAINING STATOR 
RELAXATION TIMES WITH TOTAL TEMPERATURE = 2200 °K 







AND TOTAL 


PRESSURE = 


15 ATM 




A/A* 


T 


P 


x Eqn . (22) 


T Graphs 


X Eqn. (23) 


1 


1947 


11.06 


13.3 


63,3 


61.1 


1.14 


1740 


6.38 


33.7 


172.4 


207.5 


1.27 


1656 


5.05 


50.6 


238 


353 


1. 37 


1604 


4.34 


65.8 


** 


498 


1.48 


1563 


3.85 


81.3 


** 


658 


1.51 


1544 


3.60 


90.8 


** 


759 


1.55 


1531 


3.44 


98 


** 


837 



** Unable to extrapolate the graph further. 
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TABLE V 



STATOR RELAXATION TIMES IN MICRO-SECONDS 
OF NITROGEN MOLECULES UNDER VARIED TOTAL PRESSURES 
WITH TOTAL TEMPERATURE = 2200°K USING EQUATION (23) 



A/A* 




Pressure, Atm 




20 15 


15 


10 


5 








** 






1 




61.1 82.9 


10. 8 


12 4 


249 


1.14 




207.5 277 


21 


415 


830 


1.27 




353 471 


28 


706 


1412 


1.37 




498 663 


33 


996 


1992 


1.46 




658 876 


39 


1316 


2631 


1.51 




759 1011 


42 


1577 


30 35 


1.55 




837 -1116 


44 


16 73 


3347 




** Collision Partner is H„0 in this column vice N. 


2 




as 


in the other columns. 










TABLE VI 






COMPOSITE 


RELAXATION TIMES CALCULATED USING EQUATION 


(30a) 






WITH TOTAL PRESSURE = 15 ATM 






A/A* 


n 2 -n 2 


h 2 o-n 2 


Composite 




1 


82.9 


10.8 


44.2 






1.14 


277 


21 


106.4 






1.27 


471 


27.9 


152.5 






1. 37 


663 


33.4 


190.6 






1.46 


876 


38.7 


228 






1.51 


1011 


42 


250.5 






1.55 


1116 


44.3 


267 
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TABLE VII 



RATIO OF ACTUAL VIBRATIONAL ENERGY 
TO EQUILIBRIUM VIBRATIONAL ENERGY AT THE THROAT 
UNDER VARIOUS TOTAL PRESSURES AT TOTAL TEMPERATURE = 2200 °K 



Total 

Pressure, ATM 


E* 

V 


A/A* 


E 

V 


Total 

Pressure 


A/A* 


E 

V 


E* 

V th 


E* 

v th 


E* 

V th 


20 


1 


1 


1 


10 


1 


1 




. 8388 


1.14 


.9919 




1.14 


.9958 




.7733 


1.27 


.9779 




1.27 


.9886 




.7331 


1.37 


.9634 




1. 37 


.9831 




.7006 


1. 46 


.9553 




1.46 


.9786 




.6856 


1.51 


.94 85 




1.51 


.9749 




.6755 


1.55 


.9423 




1.55 


.9715 


15 




1 


1 


• 5 


1 


1 






1.14 


.9938 




1.14 


.9979 






1.27 


.9832 




1.27 


.9942 






1.37 


.9751 




1.37 


.9914 






1.46 


.9687 




1. 4 6 


.9891 






1.51 


.9633 




1.51. 


.9871 






1.55 


.9584 




1?55 


.9854 
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TABLE VIII 



RATIO OF ACTUAL VIBRATIONAL ENERGY TO EQUILIBRIUM VIT RATIONAL 
ENERGY AT THE THROAT WITH TOTAL PRESSURE = 15 ATM AND TOTAL 
TEMPERATURE = 2200°K FOR VARIOUS MIDSPAN STREAMLINE LENGTHS 



E* 

V 


A/A* 


E 

V 






E 

V 


E 

V 


E* 

v th 


E* 

V th 






E* 

V th 


E* 

v th 






+ 






tt 


ttt 


1 


1 


1 






1 


1 


.8388 


1.14 


.9938 






.9979 


.9946 


.7733 


1.27 


.9832 






.9942 


.9846 


.7731 


1. 37 


.9751 






.9914 


.9758 


.7706 


1.46 


.9687 






.9891 


.9680 


.6856 


1.51 


.9633 






.9871 


.9609 


.6755 


1.55 


.9584 






.9854 


.9544 


t : 


Midspan streamline 

N . 


length 


is 


6"; 


collision 


Partner 


tt : 


Midspan streamline 

n 2 . 


length 


is 


2"; 


collision 


Partner 


ttt : 


Midspan streamline 


length 


is 


2"; 


collision 


Partners 



N 2 , C0 2 and H 2 0. 



TABLE IX 

CHEMICAL EFFICIENCY CALCULATED USING EQUATION (34) 
WITH TOTAL TEMPERATURE = 2200 °K 



Turbine Exit 
Total Temperature 


n c 


1500 


98.92 


1600 


99.00 


1700 


99.15 


1800 


99.15 


1900 


99.50 
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TABLE X 



ELEMENTARY CHEMICAL REACTIONS FROM NEWHALL f 10 



1. 


h 2 o 


+ 


M 


2. 


H 2 


+ 


M 


3. 


NO 


+ 


M 


4. 


N 2 


+ 


M 


5. 


no 2 


+ 


M 


6. 


°2 


+ 


M 


7. 


OH 


+ 


H 


8. 


OH 


+ 


0 


9. 


OH 


+ 


H 2 


10. 


OH 


+ 


OH 


11. 


CO 


+ 


OH 


12. 


NO 


+ 


0 


13. 


NO 


+ 


N 


14. 


NO 


+ 


°2 



<- 


OH 


+ 


H 


+ 


M 


-+ 


H 


+ 


H 


+ 


M 


<- 


N . 


+ 


0 


+ 


M 


-> 


N 


+ 


N 


+ 


M 


->• 

4r 


NO 


+ 


0 


+ 


M 


-»* 


0 


+ 


0 


+ 


M 


-*■ 

<- 


H 2 


+ 


0 






•* 

<- 


°2 


+ 


H 






-»• 

<- 


h 2 o 


+ 


H 






<- 


h 2 o 


+ 


0 






-* 

+- 


C°2 


+ 


H 






-»■ 

<- 


o 

N> 


+ 


N 






-+ 

<- 


N 2 


+ 


0 






->■ 

4- 


no 2 


+ 


0 
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TABLE XI 



SELECTED RATE CONSTANTS 
FOR 

THE ELEMENTARY CHEMICAL REACTIONS FROM NEWHALL [10] 



Reaction 



Forward Rate Constant 



Backward Rite 
Constant 



1 

2 

3 



5.4 x 10 17 exp ( - 12 3 , 6 0 0 /RT ) 
3 . Ixl0 15 exp (-110 , 000/RT) 



7x10 10 (150, OOP) 2 



,1/2 



4.2x10 



(RT) 



12 



J72 



( 22^900.) ex? { _ 2 2 4 , 9 0 0 /RT ) 



1.5 x 10 



16 



7.0 x 10 17 / T 
14 



exp (-150 , 000/RT) 9.0 x 10 



6.1 x 10 



14 



5 


2 1 

-• 4 |- 1Q exp (-74 , 000/RT) 


2.0 x 10 16 


6 


6x10 13 1;L ^ 000 exp (-11 , 800/RT) 


1.0 x lO 3- ^ 


7 


1. 4x10 ^exp (-6 , 000/RT) 


3 . 3x10 12 exp (-8 , 000/RT) 


8 


5.5xl0^ 2 exp(-l, 000/RT) 


7 . 2x10 3 ^exp (-16 , 900/RT) 


9 


6 .2xl0 13 exp (-6 , 000/RT) 


3. 2x10 14 exp (-21, 100/RT) 


10 


7.7x10 3-2 exp ( -1 , 000/RT) 


8. 3x10 13 exp( -18, 100/RT) 


11 


7.1x10 12 exp (-7 , 700 /RT) 


4 . 7x10 14 exp (-27, 250/RT) 


12 


3.2xl0 9 exp (-39 , 100/RT) 


13 . 3x10 9 exp (-7 , 0 80/RT) 


13 


1 . 55x10 13 


7x10 13 exp (-75 , 500/RT) 


14 


1.8xl0 10 T 1/2 exp(-47, 000/RT) 5.8xl0 10 T 1/2 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 



TABLE XII 



RATE CONSTANTS AT 1947°K 



Forward Constant Backward Constant 



7.1979 x 10 3 


1.5 x 10 16 


1.3896 x 10 3 


3.5933 x 10 14 


3.4565 x 10" 5 


9.0 x 10 15 


6.1003 x 10 ~ 10 


6.1 x 10 14 


1. 36733 x 10 10 


2.0 x 10 16 


1.037 x 10 2 


14 

1.0 x 10 


2.9688 x 10 11 


4.17295 x 10 11 


4.247 x 10 13 


9.1234 x 10 12 


1.3147 x 10 13 


1.3693 x 10 12 


5.9461 x 10 12 


7.7125 x 10 11 


9.702 x 10 11 


4.1025 x 10 11 


2.542 x 10 8 


4.1536 x 10 12 


13 

1.55 x 10 


2.3418 x 10 5 


4.205 x 10 6 


2.559 x 10 12 
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TABLE XIII 



EQUILIBRIUM MOLE FRACTIONS WITH TOTAL TEMPERATURE = 2200 °K 
AND TOTAL PRESSURE = 15 ATM FROM BANES, ET AL. , [3] 





Throat 


Exit 


T 


194 7°K 


1531°K 


P 


8 . 3 atm 


2.58 atm 


co 2 


1.2919xl0 -1 


1. 3073xl0 -1 


1 O 
(N 

a 


1.2992xl0 -1 


1.3048X10 -1 


°2 


6.7092xl0 -4 


5. 3073xl0~ 5 


H 2 


3. 4679xl0 -4 


3.6215X10 -5 


N 2 


7.2896xl0 _1 


7 . 29 81xl0 _1 


N 


lxlO -12 


lxlO" 12 


NO 


3.876 8xl0~ 4 


2.4451xl0 -5 


OH 


2 . 5755x10 ~ 4 


1.276 8xl0 -5 


CO 


1. 409 8x10 ~ 3 


1.006 8x10 ~ 4 


no 2 


1.1142xl0 -7 


2 . 80 80x10 -9 


H 


7.9334xl0~ 6 


1.2395xl0 -7 


0 


4.4208xl0 -6 


3 . 8529x10 


M 


28. 876 


28.903 


P 


1.499xl0 -3 


5.93 8xl0~ 4 
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TABLE XIV 



EQUILIBRIUM MOLE FRACTIONS WITH TOTAL TEMPERATURE = 2200 °K 
AND TOTAL PRESSURE = 5 ATM FROM BANES, ET AL . , [3] 





Throat 


Exit 


T 


194 7 °K 


1531°K 


P 


2.7655 atm 


0.86 atm 


n 

o 

to 


1.2845xl0 _1 


1.306 8x10 -1 


h 2 o 


1. 296 3xl0 _1 


1 . 3045xl0 _1 


°2 


1.0156xl0 -3 


7.6494xl0 -5 


H 2 


5 . 0574xl0~ 4 


5.062 8xl0~ 5 


N 2 


7 . 2 857xl0 _1 


7 . 2979xl0 _1 


N 


lxlO -12 


1x10 ” 12 


NO 


4 . 7613x10 ~ 4 


2.9394xl0 -5 


OH 


3.8269xl0 -4 


1 . 821 lxl0 _5 


CO 


2.1936X10 -3 


1.7399xl0 -4 


no 2 


9 . 4140xl0 -8 


2.4052xl0 -9 


H 


1 . 7402x10 _5 


2.4370xl0 -7 


0 


9 . 8841xl0 -6 


7 . 6 885x10 ” 3 


M 


28.863 


28.902 


F 


4.997xl0 -4 


1.98xl0 -4 
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TABLE XV 



REACTIONS RATES OF ELEMENTARY EQUATIONS AT 1947°K 



Total Pressure, Atm 


15 


5 




Forward 


Backward 


Forward 


Backward 


. 1 


1. 839xl0~ 6 


3.132xl0~ 6 


2 . 037xl0~ 7 


3. 775xl0“ 7 


2 


9 . 48xl0~ 10 


2.3127x10 


1 . 5 34xl0~ 10 


4 . 114xl0 _1 ° 


3 • 


2 . 636xl0~ 17 


4.066xl0“ 16 


_ I O 

3.593x10 


3. 36xl0 -16 


4 


8 . 748xl0~ 19 


1.02xl0~ 37 


-20 

9.703x10 


2.305xl0~ 24 


5 


2 . 99 7xl0~ 6 


3.503xl0~ 6 


2 . 81xl0~ 7 


3. 56xl0 -7 


6 


1 . 369xl0 _1 ° 


1. 997xl0~ 10 


2.30xl0 -11 


3.69xl0 -11 


7 


1. 6 37xl0~ 6 


1. 72 7x10 ~ 6 


5.92xl0~ 7 


6 . 25x10 ~ 7 


8 


1. 305xl0~ 4 


1. 3 1x10 ~ 4 


4 . 81x10 ~ 5 


4 . 83x10 ~ 5 


9 


3 . 169xl0~ 3 


3 . 809xl0~ 3 


7 . 62x10 ~ 4 


9 . 26xl0~ 4 


10 


1.064xl0“ 3 


1 . 1955xl0~ 3 


2 . 6 lxlO -4 


2.2 8xl0~ 3 


11 


1.018X10 -3 


1.1348X10 -3 


2.44xl0~ 4 


2 . 75xl0~ 4 


12 


1.176xl0~ 9 


-12 

7.52x10 


3 . 58xl0 -10 


1.26xl0 -12 


13 


1.620X10 -11 


2 . 04xl0~ 9 


-12 

2.21x10 


5 . 05xl0~ 10 


14 


2 . 952xl0~ 9 


3 . 40x10 ~ 9 


6.09xl0 -10 


7 . 14xl0~ 10 
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TABLE XVI 



EQUATIONS FOR TOTAL REACTION RATE OF EACH SPECIE 



Rc0 2 




r llf _r llb 


^2° 


= 


r lb" r lf +r 9f _r 9b +r 10f" r llb 


R °2 


= 


r 6b" r 6f +r 8f" r 8b +r l2f“ r 12b +r 14b" r 14f 




= 


r 2b -r 2f +r 7 £ -r 7b +r 9b -r 9 £ 


""2 


= 


r 4b _r 4f +r 13f ~ r l3b 


^0 


= 


r 3b -r 3f +r 5f -r 5b +r i2b -r 12f +r 13b _r 13f +r 14b -r 14f 


r oh 


= 


r lf -r lb +r 7b _r 7f +r 8b“ r 8f +r 9b“ r 9f +2r 10b" 2r 10f +r llb" r llf 


R co 


= 


r llb " r llf 


**°2 


= 


r 5b" r 5f +r 14f- r 14b 


*H 


= 


r lf- r lb +2r 2f- 2r 2b +r 7b- r 7f +r 8f- r 8b +r 9f" r 9b +r llf _r llb 


R o 


= 


r 3f" r 3b +r 5f" R 5b +2r 6f" 2r 6b +r 7f" r 7b +r 8b“ r 8f +r 10f" r 10b 






+r 12b _r 12f +r 13f" r 13b +r 14f“ r 14b 



55 



i 










TABLE XVII 



TOTAL REACTION RATE, R, AT THE THROAT 
WITH STATIC TEMPERATURE = 1947°K 



Species 


Total 


Pressure 


15 


5 


C°2 


- 1.16 8xl0 -4 


- 3. 13xl0 -5 


H 2° 


- 7 . 7xl0~ 4 


- 2 . 185xl0~ 3 


°2 


- 4.98xl0 -7 


- 1.695xlO" 7 


H 2 


6. 4xl0" 4 


1.63X10 -4 


N 2 


- 2.02xl0 _9 


- 5 . 03xl0~ 10 


NO 


- 5.047X10 -7 


- 7.445X10 -8 


OH 


1.019X10 -3 


4 . 24xl0 _3 


CO 


1. 168xl0 _4 


3 . 13xl0~ 5 


no 2 


5 . 056x10 ~ 7 


7 . 46x10 ~ 8 • 


H 


- 7.59xl0 -4 


- 1.95X10 -4 


0 


- 1. 316xl0~ 4 


- 2.022xl0 _3 
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TABLE XVIII 



CHARACTERISTIC REACTION TIMES IN THE THROAT AREA 
IN ORDER TO REACH EXIT EQUILIBRIUM CONDITIONS 
FOR TOTAL TEMPERATURE = 2200 °K 

Total Pressure 
15 5 

Species Time 



C0 2 


+ 


6.3x10 4 
-5 


+ 


1.14 

-3 


h 2 o 


+ 


2.95x10 3 




5.1x10 


°2 




6 . 4 4xlO*" 2 


- 


9 . 59x10 1 


H 2 


+ 


2.52xl0~ 5 


+ 


4.84xl0 -2 


N 2 


+ 


4.32 

-2 


+ 


8 . 10x10 3 
2 


NO 




3.74x10 

-5 




1.04x10 

-3 


OH 


+ 


1.25x10 3 
-4 


+ 


1.49x10 

1 


CO 


+ 


6.26x10 

-5 


+ 


2.00x10 

-1 


no 2 


+ 


1.12x10 3 
-7 


+ 


3.62x10 

-2 


H 




5.34x10 




2.56x10 


0 




1. 73xl0“ 6 




1.40xl0~ 3 



Transit Time 5.7x10 
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TABLE XIX 



FREEZING CLASSIFICATION OF CHEMICAL SPECIES 



Total 

Pressure 



Frozen 


Unknown 


CN 

O 


H 2° 


N 2 


co 2 


NO 


H 2 


H 


OH 


0 


CO 




no 2 


h 2 o 


None 


C°2 




°2 




H 2 




N 2 




NO 




OH 




CO 




CN 

O 

2 




H 




0 





Equi librium 



15 



None 



None 
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APPENDIX B 




FIGURE I TURBINE SCHEMATIC AND FLOW STATION 
IDENTIFICATION FOR THE VIBRATIONAL 
FREEZING ANALYSIS 
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FIGURE 2 ENERGY DISTRIBUTION IN THE TURBINE 
FOR THE VIBRATIONAL FREEZING ANALYSIS 
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Midspon Streamline 




FIGURE 3 STATOR PASSAGE 
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CHANGE IN TOTAL TEMPERATURE ACROSS THE 
TURBINE, T tjj - Tf 2 , °K 



FIGURE 4 PERCENTAGE PERFORMANCE LOST DUE 
TO VIBRATIONAL FREEZING 
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FIGURE 5 LOCAL VIBRATIONAL RELAXATION TIMES FOR 
THE STATOR PASSAGE IN FIGURE 3 
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AREA RATIO, A / A* 



FIGURE 6 RATIO OF LOCAL VIBRATIONAL ENERGY TO 
EQUILIBRIUM THROAT VIBRATIONAL ENERGY 
AS A FUNCTION OF AREA RATIO FOR 
PURE NITROGEN 



64 




FIGURE 7 INFLUENCE OF STATOR SIZE AND COLLISION 
PARTNER ON THE RATIO OF LOCAL VIBRATIONAL 
ENERGY TO EQUILIBRIUM THROAT VIBRATIONAL 
ENERGY AS A FUNCTION OF AREA RATIO 
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COMPLETE EQUILIBRIUM 



ENTHALPY 




ENTROPY 




FIGURE 8 EXPANSION PROCESSES USED IN 
ANALYSIS OF CHEMICAL FREEZING 
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CHANGE IN TOTAL TEMPERATURE ACROSS THE TURBINE 




FIGURE 9 PERCENTAGE PERFORMANCE LOST AS A 
RESULT OF CHEMICAL FREEZING 
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APPENDIX C 



An approximate solution to equation (21) may be obtained 
using the finite difference method. Equation (21) is put in 
the form as in Hurle, et al. , [9], 

[1 - (Ax/2X )] 

E ^ T v^n+l - [1 + (Ax/2X n ^) ] E ^ T v^n 

Ax/2 E(T *n+l , E (T) n 

[l + (Ax/2 W l Vl A n 

where n and n + 1 refer to two adjacent locations along the 
nozzle axis separated by distance Ax, and X and T are evalu- 
ated for equilibrium flow. X is the characteristic relaxation 
length and X = Ux, where U is the local gas velocity and t is 
relaxation time. 
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APPENDIX D 



This is the calculation of the conversion constant for 
use in the rate equations so that mole fractions may be 
utilized directly. 

The units of the rate constants in Table XI vary with 
the type of elementary equation; bi-mole cular or tri- 
molecular. The normal units for concentrations are moles 
- cm and the units of X are moles of specie i per mole of 
mixture. From the ideal gas law 

P V = n R T 

Multiplying by moles of specie i over moles of specie i 



P V X. = n. R T 
1 i 



or 



n i ' 



P X. 
1 



V R T 

_3 

in units of moles - cm ; the desxred units, 
Using the gas constant equal to 



R = 82.06 



cm atm 



mole °K 

a conversion factor for use at a specified temperature and 
mixture pressure can be obtained. The mole fractions are each 
multiplied by this factor and are then ready for use in the 
specie rate equations. 
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